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The radiatively induced neutrino mass matrix is investigated on the basis of an SU(5) SUSY 
model. In order to evade the proton decay, an ansatz based on a discrete symmetry Z2 is assumed: 
although, at the unification scale, we have two types of superfields 9l(±) = ^l(±) + Wl(±) , which are 
transformed as ^l(±) — * i^£(±) under the discrete symmetry Z2, the particles are decoupled 

after the SU(5) symmetry is broken, so that our quarks and leptons belong to The _R-parity- 

violating terms for our quarks and leptons ^l(-) are basically forbidden under the symmetry Z2. 
However, we assume that mixings between members of ^l(+) an d those of ^l(-) are m part caused 
after SU(5) is broken. As a result, the J?-parity-violating interactions are in part allowed, so that 
the neutrino masses are radiatively generated, while the proton decay due to the J?-parity violating 
terms is still forbidden because the term d^d^u^ has z = — 1. 

PACS number(s): 11.30.Er; 12.60. Jv; 14.60.Pq; 11.30.Hv 



I. INTRODUCTION 

The origin of the neutrino mass generation is still a 
mysterious problem in the unified understanding of the 
quarks and leptons. The Zee model [1] is one of several 
promising models, because it has only 3 free parame- 
ters and it can naturally lead to a large neutrino mixing 
[2], especially, to a bimaximal mixing [3]. However, the 
original Zee model is not on the framework of a grand 
unification theory (GUT). The most attractive idea [4] 
to embed the Zee model into GUTs is to identify the 
Zee scalar h + as the slepton e R in an i?-parity-violating 
supersymmetric (SUSY) model. However, usually, it is 
accepted that SUSY models with i?-parity violation are 
incompatible with a GUT scenario, because the i?-parity- 
violating interactions induce proton decay [5]. 

In the present paper, in order to suppress this kind 
of proton decay, a discrete symmetry Z2 is introduced. 
The essential idea is as follows: At the unification scale 
fi = Mx, we have two types of superfields ^l(±) — 
5l(±) + 10l(±)> which are transformed with ±1 under the 
discrete symmetry Z2 (we will call it "Z2-parity" here- 
after). We consider that the particles are decou- 
pled after the SU(5) symmetry is broken (but Z 2 is still 
unbroken), so that our quarks and leptons 5l + 10l be- 
long to \& The i?-parity violating terms are given 
by the combinations 5( + )5( + )10( + ), 5(_)5(_)10( + ) and 
5(_|_)5(_)10(_), so that they basically do not contribute 
to the quarks and leptons with Z2-parity z = — 1, be- 
cause of the Z2 symmetry. However, we assume that 



mixings of the members of with those of ^l(-) 

are caused in part after SU(5) is broken. As a result, 
the i?-parity- violating interactions are in part allowed, so 
that the neutrino masses are radiatively generated, while 
the proton decay due to the i?-parity-violating terms is 
still forbidden because the term d c R d c R u R is still exactly 
forbidden below /1 = Mx in the present scheme. The 
details will be discussed in the next section. 

The purpose of the present paper is to investigate the 
possible forms of the radiatively induced neutrino mass 
matrix under the Z2 symmetry. In Sec. Ill, we will 
give them, including a numerical study. In Sec. IV, we 
will give a comment on the Higgs scalars in the present 
scheme. Finally, Sec. V is devoted to our conclusion. 

II. Z 2 SYMMETRY AND THE PROTON DECAY 

We identify the Zee scalar h + as the slepton e^, which 

is a member of SU(5) 10-plet sfermions ipio. Then, the 
Zee interactions correspond to the following i?-parity- 
violating interactions 

= -j=\ ijk {e aM {d B )«{d c R f 3 (rf R )l 

~[(e c L )i^L) 3 - {vl)i{e L )j]{e R )k 
-[(Si)i(^-(rffl)?(e i ) j ](S i ) fca 

+ [PL)i(*k)f ~ , (2.1) 
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where ip c = Cip and the indices (A,B,---) 
and (a,/3, ■■■) are family-, SU(5)ge/t- and SU(3) co / our - 
indices, respectively. The coefficients are antisym- 
metric in i and j. On the other hand, in SUSY GUT 
models, if the interactions (2.1) exist, the following im- 
parity-violating interactions will also exist: 

= -^ Xi 3 k { £ "/37(^)?(^L)i( u fl)fc 

-[(ei)i(^L)j - {v C LWL) 3 }(e R ) k 

-[&)A)j-(dR)i@L)i](UL)ka 

+[(^)i(d t «) J ? -(dfl)?(^) J -](dL)fca} , (2-2) 

which contribute to the proton decay through the inter- 
mediate state cIr. Also, the term {d R )f{d c R )^ (u R )1 in 
the interactions (2.1) can contribute to the nucleon de- 
cay through the intermediate state u R . The upper limits 
of the coupling constants \j k from proton decay exper- 
iments have been investigated by Smirnov and Vissani 
[5], and the values must be highly suppressed. 

In order to forbid the contribution of the interactions 
(2.1) and (2.2) to the proton decay, we must consider that 
in the i?-parity-violating interactions 5 x 5 x 10, the term 
d c R d R u c R is exactly forbidden, while the terms vl&l& r 
and/or vi,d c R dh are in part allowed. 

For such purpose, we introduce a discrete symmetry 
Z2, which exactly holds at every energy scale. At the 
unification scale \i — Mx, we have two types of super- 
fields *l(±) = 5i,(±) + 10l(±), which are transformed 
as *l(±) — * ±^l(±) under the discrete symmetry Z 2 . 
We consider that the particles are basically de- 

coupled after the SU(5) symmetry is broken, so that our 
quarks and leptons (and their SUSY partners) 5l + 10 l 
are regarded as ^l(-) = ^l(-) + 10 l(-)- The i?-parity- 
violating terms for quarks and leptons (and their SUSY 
partners) are basically forbidden under the symmetry 
Z 2 below /j, = Mx, because the terms are composed of 
5 i( _)5 i( _ ) 10 L( _ ) . 

However, if wc assume that mixings between the mem- 
bers of an d those of ^l(-) m part take place 
after SU(5) is broken, i?-parity-violating interactions 
^ l<-)^ <-) become available at the low energy 
/i = mz, too. For example, we assume a mixing 

(2, l) Li = (2, l) L( _ }i cos Of + (2, l) L(+)l sin 0f , (2.3) 

between the (2,1) components of SU(2)xSU(3) for the 
i-th family. (Hereafter, we will refer to the mixing (2.3) 
as a mixing of type Aj.) Then, the i?-parity- violating 
interactions 

sin 8 <f\ijkVLid c Rj d L k , sin 8 >f\ijkeLid c Rj u L k , 

sin 6f cos 8f \ijkVLie L je c Rk , (2.4) 

become available from the interactions 



Ay7 s 5( + )j5(_) j 10(_) fe , (2.5) 

above the unification scale /j, = Mx- Also, we can con- 
sider a mixing 

(2, 3) Lfe = (2, 3) L( _ )fc cos 6* + (2, 3) L(+)fc sin fljf , (2.6) 

between the (2,3) components of SU(2)xSU(3) for the 
k-th family. (Hereafter, we will refer to the mixing (2.6) 
as a Bfc-type mixing.) Then, the i?-parity- violating in- 
teractions 

sin(?f X' ijk v Li d R jd L k and sinflf X' ijk e L id c Rj ULk (2.7) 

become available from the interactions 

^fc5(_) i 5(_) j 10( + ) fe . (2.8) 

On the other hand, note that the interaction 

d R d R u R (2.9) 

is exactly forbidden, independently of whether the mix- 
ings (2.3) and (2.6) occur or not, because those interac- 
tions have the Z 2 parity z = — 1. Therefore, the proton 
decay due to the i?-parity-violating terms is exactly for- 
bidden because of the absence of the term d c R d c R u R . On 
the other hand, the neutrino masses are radiatively gener- 
ated through the interactions v^d c R dh and vlGlg c r with 
z = +1. The possible forms of the radiative neutrino 
mass matrix will be discussed in the next section. 

At present, we do not know a reasonable mechanism 
not only for such a mixing, but also for the decoupling of 
*L(+). In order to make */,(+) = 5 L(+) + 10 L(+) heavy, 
the SU(2)l symmetry must be broken, but, of course, 
we cannot consider a scenario in which SU(2) is broken 
just after SU(5) is broken. In the present paper, we give 
only a phenomenological selection rule: if the superfield 
W can make a five-body SU(5) singlet operator ^ \I> \I/ \I> 
with the Z 2 parity z = +1, then the superfield \f can 
be decoupled below n = Mx- Obviously, according to 
this selection rule, the superfield can be decoupled 

below /i = Mx- Similarly, the superfield 10u + \ is decou- 
pled below (i = M X - However, note that those operators 
in the SU(5) singlets are symbolically expressed in terms 
of SU(2)xSU(3) components as follows: 

(5 i(+) ) 5 = [(2, l) i(+) ] 2 x [(1, 3) L(+) ] 3 , (2.10) 
(lO i(+) ) 5 = (1,1) M+ ) x [(l,3) L(+) ] 2 x [(2,3) i(+) ] 2 

+ (l,3) i(+) x [(2,3) L(+) ]\ (2.11) 

and that even if the interchanges (2, 1)l(+)% <-> 
(2,l) L( _ )i , and/or (2,3) L(+)fe <-> (2,3) L( _ )fc , are caused, 
the composite operators 

(5 l ) 5 = [(2,1)l(-)] 2 x[(M) l(+) ] 3 , (2.12) 
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(10 L ) 5 = (l,l) L(+) x [(l,3) i(+) ] 2 x [(2,3) iM ] 2 

+ (l,3) L(+) x [(2,3) i( _)] 4 (2.13) 

still have z = +1. Such interchanges are possible only 
for the components (2, 1)l and (2,3)l- As a result, only 
the combination 

5 L + 10 L = [(2,l) i(+) + (l,3) i( _)] 
+ [(l,l) i( _ ) + (2,3) i( _ ) + (l,3) L( _ ) ] , (2.14) 

for the i-th family and/or 

5 L + 10 L = [(2, I)l(-) + (!,%(-)] 
+[(l,l) L( _ ) + (2,3) L(+) + (l,3) L( _ ) ] , (2.15) 

for the k-th family survive below /x = Mx as the quarks 
and leptons (and their SUSY partners). 

Of course, the above selection rule cannot be justi- 
fied within the framework of the minimal SUSY standard 
model. At present, this is only an ansatz to select which 
components of SU(2)xSU(3) can be interchanged. 



III. RADIATIVE NEUTRINO MASSES 

In a SUSY GUT scenario, there are many origins of 
the neutrino mass generations. For example, the sneu- 
trinos can have vacuum expectation values (VEVs), 
and the neutrinos VLi acquire their masses thereby (for 
example, see Ref. [6]). Although we cannot rule out a 
possibility that the observed neutrino masses can be un- 
derstood from such compound origins, we do not take 
such a point of view in the present paper, because the ob- 
served neutrino masses and mixings appear to be rather 
simple and characteristic. We simply assume that the ra- 
diative masses are only dominated even if there are other 
origins of the neutrino mass generations. 

In the present scenario, the origins of the radiatively 
induced neutrino masses are two: one is induced by the 
i?-parity-violating interactions v^d c R d^ and v^d R d^ the 
other one is induced by vi,ejji R and vj,ei J e R . Note that 
there is no Zee-type diagrams due to H^-e R mixing in 
this scheme. 

First, we discuss the down-quark loop contributions. 
For simplicity, we assume that the masses Mn and Mm 
of the squarks dn and dm are approximately constant, 
independently of the flavours, although we consider the 
flavour-dependent structure for the mass terms d} L M^dm 
Then, the radiatively induced neutrino mass matrix due 
to the A-type mixing is given by 

{M u )ij = m \ lkm \ jln (Ml) kn (M^) lm + (i <-> j) , (3.1) 

where sine-factors have been drooped for simplicity, M v 
is defined by VlM v v c l , and the coupling constants Xij k 
are redefined by 



VLidRjd} Lk + v Li d Rj d° Lk - (v L <-> d° R ) . (3.2) 



Here, we have changed the definition of \j k from that 
in (2.1) as Xij k — * X*j k for the convenience of the ex- 
pression of M v defined by v^M v v c L . In the present pa- 
per, the unitary matrix U v used to diagonalize the Ma- 
jorana mass matrix M v is defined as UlM u U* = D v . 
Then, the so-called Maki-Nakagawa-Sakata-Pontecorvo 
[7] matrix (we will simply call it the "lepton mixing ma- 
trix") U = Umnsp is given by U = U^U V . Usually it 
is considered that the matrix form of M 2 is proportional 
to the form Md- Then, the neutrino mass matrix (3.1) 
becomes, in a more concise form: 



{M v )ij = m \ ikm \ji n (M\) kn (Ml)i 



(3.3) 



where we have redefined the common factor m from that 
in (3.1). Of course, exactly speaking, for the Aj mixings 
with sin ^-factors, we should read the expression (3.3) 
as 

(M v )ij =m stsfXi km \jin{Ml) kn {Ml)i m , (3.4) 

where sf = sin Of defined in (2.3). When we consider 
the Bfe mixings, we read the expression (3.3) as 

{M u ) l j =mos^s^X' ikm X' jln (Ml) kn (Ml)i m , (3.5) 

where s k = sm8j^ defined in (2.6). For mixed-type mix- 
ings of A, and Bfe, we read (3.3) as 

{M v )ij = m (sfc^X ikm + cfs^X' ikm ) 

(sfc*X jln + cf S y jln )(M\) kn (Ml) lm , (3.6) 

where cf — cos Of and c k — cos 0^ . 

The contributions from the charged lepton loops arc 
essentially the same as (3.3), except for the absence of 
the B-type mixing and the replacement Md — > Mj . For 
simplicity, we will continue the investigation for the case 
of the down-quark loop contributions. 

For the phenomenological study of the mass matrix 
(3.3), it is convenient to take the basis on which the 
down-quark mass matrix Md is diagonal: 



U?M d U R = D d = diag(r<, ml m a 3 ) 



(3.7) 



We consider that, on the basis with Md = Dd, the 
charged lepton mass matrix M e is also approximately di- 
agonal, Vjf M e U R ~ D e = diag(raf , m|, m§), so that the 
unitary matrix U v approximately gives the lepton mixing 
matrix U = U^U V . Then, we can express (3.3) as 

(M„)n - (m d 3 ) 2 (X 133 ) 2 + K) 2 (A 122 ) 2 + 2mfm 2 ! A 123 A 132 
(M„) 2 2 = {™ d 3 ) 2 {X 233 ) 2 + K) 2 (A 211 ) 2 + 2mfmfA 2 i 3 A 23 i 
(M„)33 = (m 2 ) 2 (A 322 ) 2 + K) 2 (A 311 ) 2 + 2mXA 312 A 321 
(M„)i2 = (m3) 2 Ai 33 A 23 3 + m3m 2 Ai 23 A 232 

+mimf Ai3iA 2 i 3 + TO 2 TO^Al 2 l A 2 1 2 , 
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(M„)i 3 = (m^) 2 Ai22A 3 22 + mfm^Ai32A 3 23 

+m3mjAi3iA 3 i3 + m2m^Ai2iA 3 i2 , 

(M„) 23 = (mi) 2 A 2 iiA 3 ii + TO 3 m2A 2 32A 3 2 3 

+mfmf A 23 iA 3 i 3 + m^nii A212 A 3 2i , 



where, for simplicity, we have dropped the common factor 
tuq. In order to give the best-fit values for the observed 
neutrino data [8,9] 



R: 



Am 2 2 



5.0 x l(T 5 cV 2 
2.5 x 10- 3 eV 2 



= 2.0 x 10" 



(3.9) 



M v ~ (,f ) 2 K) 



£(^3 
14 



r-2) 
er 3 



e(r 3 +r 2 ) l + er 3 (l + er 2 ) 2 



(3.14) 



(3-8) wnere e — s 3 X/sfX' «l,r 2 = mf/mf and r 3 = TO3/7 



and we have assumed ~ 1. This mass matrix can 
give a reasonable value of R together with a nearly bi- 
maximal mixing. For example, for the parameter value 
e = 0.000374, we obtain the following numerical results: 



m\ = 0.0822to , m u 2 = -0.3276m , 



m u 3 = 2.2604m , 
(3.15) 



sin 2 26 solar = sin 2 26 l2 = 0.76 , 

[tan 2 9 solar = 0.34] , (3.10) 

sin 2 20 atm = sin 2 20 23 = l.O, (3.11) 

we must seek a parameter set that gives (M 1/ ) 2 2 — {M v ) 33 
and (M v )i2 ~ (Af 1/ ) 13 for the expression (3.8). When we 
consider the A r type mixings, we obtain [10] 

(e 2 e e \ 
ell), (3.12) 
e 1 l) 

for e = sf js 2 = st/sf and X2j 3 /X 3 j2 — m 2 /TO 3 . Gen- 
erally, when we consider only Aj-type mixings, the so- 
lutions are highly dependent on the fine-tuning among 
the coefficients Xijk, and, besides, since the mass ma- 
trix is too near to a rank-1 matrix, it is difficult to give 
the solution a small but sizeable value of R (it leads to 
an extremely small value of R). Even if we take the 
charged lepton loop contributions into consideration, the 
situation is not improved unless we assume the special 
parametrization for Xijk (A2j 3 /A 3 j2 — m^/m^ and so on). 

Next, we consider the case of the Bfc-type mixings. 
When we consider a Bi mixing, we obtain a simple mass 
matrix form 

/o \ 

M v = (a? f{mif (A 211 ) 2 X' 21} X' 31 , , (3.13) 
\ A 211 A 311 (A 311 ) J 

because the case gives the relations A- j2 sf = A-^-gsf = 0. 
It is natural to consider that X' 211 ~ A 311 since those 
come from the same interactions (2.8) [not from (2.5)]. 
Therefore, the mass matrix (3.13) can give a maximal 
mixing between and v T . A small additional term will 
reasonably give a bimaximal mixing. 

For example, we consider a mixed- type, A 3 and Bi, 
mixing. In this case, the X^k in (3.3) must be re- 
placed according to Table I. It is natural to consider that 
X^k — const = A and X\- k ~ const = A'. Then, the mass 
matrix M v is reduced to the form 



/ 0.8726 -0.4822 0.0776 \ 
U= -0.3925 -0.5978 0.6990 , (3.16) 
\ 0.2907 0.6404 0.7109 / 

i.e. 

R = 0.0201, (3.17) 
sin 2 26 12 = 4C/ 1 2 1 [/ 1 2 2 = 0.708 , (3.18) 

sin 2 26> 23 = 4C/ 2 2 3 [/ 3 2 3 = 0.988 . (3.19) 

These values are in good agreement with the best fit 
values (3.9)— (3.11) for the observed neutrino data. Al- 
though it is difficult in the original Zee model to give 
a sizeable deviation of sin 2 26*i2 from 1 [11] (it must be 
sin 2 2#i2 = 1.0), the present model can give a reasonable 
deviation from sin 2 20 ±2 = 1.0. The result 

Uf 3 = 0.00602 (3.20) 

is also consistent with the present experimental upper 
limit 

\U 13 \ 2 < 0.03 , (3.21) 
from the CHOOZ collaboration [12]. 

IV. HIGGS SECTORS 

In the present model, the quark and charged lepton 
mass matrices are generated by the VEVs of the Higgs 
scalars with z = +1. Therefore, even if SU(2)l is broken 
later, the Z 2 symmetry still exactly holds. 

Let us show the mass matrices Mf for the case of the 
mixed-type mixing A 3 and B\ as an example of the ex- 
plicit forms of Mf. 

(ci bi ai \ 
C2 b 2 a 2 , (4.1) 
Cic 3 cfb 3 c^a 3 J 
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•fci 


cfc 2 


cfc 3 


h 






ai 


«2 


A3 



form, which can give satisfactory numerical results for the 
Md = [ 61 &2 63 I , (4-2) observed neutrino data. How we can justify the ansatz 

is a future task for us. 



(cfc'i cfc' 2 cfc' 3 \ 
K b' 2 b' 3 . (4.3) 
a[ a' 2 a' 3 J 

Here, Mf have been defined by J^M^r (/ = u,d,e). 
Note that the mass matrix Md has a form different from 
Mj because of the mixing factors. Usually, if we consider 
one type of Higgs scalar of SU(5) 5-plet (5-plet), it is dif- 
ficult to obtain realistic mass matrices Mf (/ = u, d, e). 
Therefore, the present model has a possibility to improve 
this problem. However, whether we can give reasonable 
mass matrix forms of Mf or not is a future task to us. 

Now, we would like to give a comment on the Higgs 
sectors. In Sec. II, we have assumed that although the 
superfield is decoupled below /j, = Mx, the com- 

ponents (2, 1)(+) of 5x,(+) can contribute to low energy 
phenomena through the mixing (2.3). If we consider the 
Higgs fields H d (+) and H d t-) with SU(5) 5-plet, and if we 
assume a situation similar to the matter fields 5l, then 
we obtain 

i? d = (2,l) (+) + (l,3) ( _) , (4-4) 

where we have assumed a perfect interchange between 
(2, 1)(_) and (2, 1)(+), not a mixing. Note that in this 
scheme, the Higgs scalar component (3, 1)(-) cannot cou- 
ple to the fermions dnu c R with z = +1 independently 
of the mixings Aj and Bk, so that the scalar (3, 1)(-) 
cannot contribute to the proton decay and it need not 
be super-heavy. (Although it couples to the fermions 
{p\dh —e c L UL), these interactions cannot contribute to 
the proton decay.) We can consider a similar mecha- 
nism for the Higgs fields H u with SU(5) 5-plet. The 
interactions of (1,3)(_) with u c L dr, and &ru c r are ab- 
sent, so that the scalar (1, 3)(_) need not be super-heavy. 
The so-called /z-terms are composed of H u i y+ )H d (+) and 
H u {-)Hd(-)- 

V. CONCLUSION 

In conclusion, we have investigated possible forms 
of radiatively induced neutrino mass matrix under an 
ansatz [(2.3) and (2.6)] within the framework of the SU(5) 
SUSY model. We have assumed two types of matter 
fields "J/ l = 5l + 10l, ^(±), which are transformed as 
^(±) — > i'J'Vi) under the Z 2 symmetry. We have as- 
sumed that the Z 2 symmetry exactly holds, even if SU(5) 
is broken. The essential ansatz is in the mixings (2.3) and 
(2.6). Although the origin of the mixings (2.3) and (2.6) 
is still an open question, if we admit this ansatz, we can 
obtain very interesting and simple neutrino mass matrix 
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TABLE I. Rule of the replacement A^fe in the mass matrix 
(3.3) for the case of A3 and Bi mixings. 
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